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<N ; Abstract 

z£\ • We study shear viscosity in weakly coupled hot cf) 4 theory using the closed 

time path formalism (CTP) of real time finite temperature field theory. We 
show that the viscosity can be obtained as the integral of a three-point func- 

<h- | tion. Although the three-point function has seven components in the CTP 

formalism, we show that the viscosity is given by a decoupled integral equa- 
tion which involves only one retarded three-point function. Non-perturbative 
corrections to the bare one loop result can be obtained by solving a Schwinger- 
Dyson type integral equation for this vertex. This integral equation represents 
the resummation of an infinite series of ladder diagrams which all contribute 
O ,' to the leading order result. It can be shown that this integral equation has 

q; exactly the same form as the Boltzmann equation. We show that the integral 

equation for the viscosity can be reexpressed by writing the vertex as a com- 
bination of polarization tensors. An expression for this polarization tensor 
can be obtained by solving another Schwinger-Dyson type integral equation. 

/\ ' This procedure results in an expression for the viscosity which represents a 

non-perturbative resummation of contributions to the viscosity which includes 
certain non-ladder graphs, as well as the usual ladders. We discuss the sig- 
nificance of this set of graphs. We show that these resummations can also be 
obtained by writing the viscosity as an integral equation involving a single 
four-point function. Finally, we show that when the viscosity is expressed 
in terms of a four-point function, it is possible to further extend the set of 
graphs included in the resummation by treating vertex and propagator cor- 
rections self-consistently. We discuss the significance of such a self-consistent 
resummation and show that the integral equations that are involved contain 
cancellations between vertex and propagator corrections. Lastly, we discuss 
the prospect of generalizing our technique to gauge theories. 
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I. INTRODUCTION 

The investigation of the transport properties of hot dense matter is a topic of great inter- 
est in the context of relativistic heavy-ion collisions and astrophysics Jl]||. The application 
of kinetic theories to the systems which are relativistic and quantum in nature involves diffi- 
culties which have been only partially overcome ||. There are two methods that are used to 
calculate transport coefficients. The first method involves solving the transport equations. 
Transport coefficients are calculated from classical transport equations using the relaxation 
time approximation [Q. The second method is to use the Kubo formulae which relate trans- 
port coefficients to the low frequency, zero momentum limits of the spectral densities of the 
appropriate composite operators []5|-§]. The Kubo formulae provide a framework for calcu- 
lating transport coefficients starting from first principles using finite temperature quantum 
field theory. 

The diagrammatic rules for computing spectral densities have been derived in both the 
real time formalism |J and the imaginary time formalism |IOj of finite temperature field 
theory. One-loop calculations of transport coefficients using Kubo formulae with effective 
propagators (including the single particle thermal lifetimes) have appeared previously ||T0| , |TT| . 
However, these calculations are incomplete, even in the weak coupling limit, because they do 
not include an infinite class of planar ladder diagrams, which contribute at the same order as 
the one-loop diagrams. The resummation of this infinite series of ladder diagrams is highly 
nontrivial and has been performed by Jeon JT^JT^] using the imaginary time formalism. 



In the imaginary formalism, this resummation is extremely complicated. In this paper we 
perform the ladder resummation in a more efficient way using the closed time path (CTP) 



framework fl3HTq[ . We give a compact derivation of the resummation of the infinite series 
of ladder diagrams that contribute to the shear viscosity. We can write the viscosity as 
an integral equation involving a single three-point function which can be obtained as the 
solution of a further decoupled linear integral equation. We show that these results are 



the same as those obtained in [10]. We will also show that the ladder resummation can 
be obtained by writing the viscosity as an integral equation involving a single four-point 
function, which is obtained, as before, from a decoupled linear integral equation. The four- 
point function is introduced because it allows us to obtain a special type of self-consistent 
resummation which includes additional non-ladder graphs. It is important to note that 
in the CTP formalism the three-point function has seven components and the four-point 
function has fifteen components, which means that the fact that these integral equations are 
decoupled is highly non-trivial. 

We also study extended resummations which include non-ladder contributions. We derive 
expressions for the three- and four-point functions in terms of a difference of self energies. 
By writing the self energy as the solution of a decoupled linear integral equation we obtain 
an extended resummation. We discuss the importance of non-ladder contributions to this 
resummation. 

Finally, we show that when the viscosity is expressed in terms of a four-point function, it is 
possible to perform a self-consistent resummation in which corrections to the vertex and the 
propagators are put on the same footing. We find there are cancellations between vertex and 
propagator corrections. For 4 theory, these cancellations do not significantly simplify the 
form of the integral equation, however, for a gauge theory with a three-point interaction, like 



scalar QED, we expect that the resulting integral equation will be considerably simplified. 
We present the motivation for this argument. In the last section we discuss our conclusions. 
This paper is organized as follows. 
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II. NOTATION 
A. The CTP formalism 



The CTP formalism of finite temperature field theory was introduced by Keldysh [|nj and 



Schwinger JTJ]. Good reviews are found in [[H],[n|[rj| . The CTP contour has two branches: 
C% runs from negative infinity to positive infinity just above the real axis, and C 2 runs back 
from positive infinity to negative infinity just below the real axis. All fields can take values 
on either branch of this contour, which results in a doubling in the number of degrees of 
freedom. We will consider 4 theory. The Lagrangian is given by, 



\ m *-\m*f-±< 



At finite temperature, when hard thermal loop corrections are included, the mass term is 
usually dropped relative to the real part of the hard thermal loop self energy. The scalar 
propagator is given by, 

D{X-Y) = -i{T c <P{X)<P(Y)) (1) 

where T c is the operator that time orders along the CTP contour. We also use the notation 
X = (t,x) and P = (po,p). The propagator has 2 2 = 4 components and can be written as 
a 2 x 2 matrix 



D 



Dxx D 12 
D21 D 22 



with 



D n (X-Y) = -i(T{<f>{X)(j>(Y))), 

D 12 {X-Y) = -i(<f>{Y)<f>(X)), 

D 21 {X-Y) = -i(<j>(X)</>{Y)), 

D 22 (X -Y) = -i{f{4>{X)4>{Y))) , (2) 

where T is the usual time ordering operator, and T is the anti-chronological time ordering 
operator. These four components satisfy, 

D n - D 12 - D 21 + D 22 = (3) 

as a consequence of the identity 6{x) + 9(— x) = 1. 

It is more useful to write the propagator in terms of the three functions 

D R = D u - D 12 , 
D A = D n - D 21 , 
D F = D lx + D 22 . (4) 

Dr and Da are the usual retarded and advanced propagators, satisfying 



D R (X - Y) - D A (X -Y) = -i([<f>{X), <f>(Y)]) , (5) 

and Dp is the symmetric combination 

D F {X-Y) = -i({<j>{X),<f>Qr)}), (6) 

which satisfies the KMS condition. In momentum space 

Dr,a(P)= 7 — t . w =^ n' 

[Po ± te) z — p z — m z 

D F (P)= (1 + 2n(p ))(D R (P) - D A (P)), (7) 

where n(p Q ) is the thermal Bose-Einstein distribution, 

n{p ) = ePpo _ - , n(-po)=-(l+n(po)), iV p = 1 + 2n{ Po ) (8) 

The propagator can be rewritten as an outer product of two component column vectors [|TH[ : 

">-° -G)(-i) + M-i)(l) + M3(l> 

Using the KMS condition (0) this expression can be rewritten, 

™-M9C^)- D H:&>)G)- 

We can extract the 1PI two-point function, or self energy, by removing external legs in 
the usual way. We find, 

n^ = iiu + n 12 
u A = n n + n 21 

u F = n n + n 22 (ii) 

where n^ and n^ are the usual retarded and advanced self energies. The symmetric com- 
ponent n^ satisfies the KMS condition, 

U F = (l + 2n(p ))(U R ~U A ) 

and the four CTP components satisfy the constraint, 

n n + n 12 + n 2i + n 22 = o 



B. Definition of Viscosity 

If a system is slightly perturbed away from equilibrium, fluctuations will occur. The 
response to these fluctuations will be characterized by transport coefficients, for instance 
shear and bulk viscosities in a system with no conserved particle number. The shear vis- 
cosity characterizes the diffusive relaxation of transverse momentum density fluctuations 
and is proportional to the two body elastic scattering mean free path. The bulk viscosity 
characterizes the departure from equilibrium during a uniform expansion and is proportional 
to the mean free path for particle number changing processes. Since the bulk viscosity in- 
volves particle number changing processes, its calculation is more difficult. However, at high 
temperature in a weakly coupled system, the shear viscosity is much larger than the bulk 
viscosity. For simplicity we will consider only the shear viscosity in this paper. 

Using standard linear response theory one may express the shear viscosity in terms of 
the stress tensor-stress tensor correlation function. One obtains the Kubo formula p|,|TU|, 



V = ^kmKmMQ) (12) 

where, 

a(Q) =JdXe lQ < x - Y \n lm (X)rc lm (Y)) (13) 

and 

n lm (X) = V z 0V m - -5 lm {V<Pf (14) 

in 4 theory. We substitute (|14"D into (|13|) and use Wick's theorem and (§). We obtain, 



a 



(Q) = 2 J dKI lm (k, k)B 12 (K)B 21 (K + Q)I lm (k, k) (15) 



This result is shown diagrammatically in Fig. [1]. The factor I[ m is associated with the 
joining of two propagators and has the form, 

Ii m (k, k) = kik m - -5 lm k 2 (16) 

where we have dropped the q dependence because we intend to take q to zero at the end of 
the calculation. In the factors I\ m no difficulties will arise from dropping the q dependence 
at the beginning of the calculation. If we attempt to work perturbatively, the lowest order 



contribution to (|15 ) is the one loop graph, as shown in Fig. [1]. For calculational purposes, 
it is easier to work with the combination of one loop diagrams shown in Fig. [2]. The first is 
a and the second we call a'. The difference of these graphs is related to a as follows. Using 
the KMS condition (|7]) in the form, 

(l + n(p ))D 12 (P)=n(p )D 21 (P) 

it is straightforward to show that, 



n(q )a = (1 + n(q ))a' (17) 

which means that we can obtain a from the difference a — a': 

a=(l + n(q ))(a-a') (18) 

It has been known for some time that there are graphs that are higher order in the loop 
expansion than the one loop graph which, nevertheless, contribute to the same order in 
perturbation theory. The purpose of this paper is to identify these graphs, and to develop 
techniques for resumming them. 

III. VISCOSITY FROM CORRECTED THREE-POINT FUNCTIONS 

A. Notation 

We define the bare three-point vertex as, 

T { Z m = lcrila5 ab 5 cb I lm (k,k) (19) 

where 1 is the two by two identity matrix and t% is the third Pauli matrix. Corrections to 
this vertex are given by the 1PI part of the three-point function, 

r2»(z, x, y) = (T c c (zKV m (x)0 a (y)> (20) 

In the real time formalism, the three-point function has 2 3 = 8 components, since each of 
the fields can take values on either branch of the contour. Only seven of these components 
are independent because of the identity 

EEE(-i) a+fe+c ' 3 rSr = o. (21) 

a=l 6=1 c=l 

The corrected three-point vertex is obtained by truncating external legs to obtain the 1PI 
part of the three-point function. These 1PI vertices obey the constraint, 

o=l 6=1 e=l 



We define the combinations, 



-plm -plm I -plm , -plm , -plm 

1 R — l 111 i L 112 + L 211 + L 212 

-plm -plm _i_ -plm , -plm _i_ -plm 

1 Ri — L 111 "r L 112 "r L 121 "+" L 122 

plm -plm I -plm , -plm , -plm 

1 Ro — L 111 "r L 121 ~r 1 2 11 ~r 1 2 21 

■plm -plm _i_ -plm _i_ -plm _i_ -plm 

1 F ~ l 111 ~r 1 12 i T 1 212 "+" L 222 

pZm -plm I pZm _i_ pirn _i_ pZm 

1 Fi — L 111 ~r 1 12 2 "I" L 211 "+" L 222 

pim pirn I -plm , -plm , -plm 

1 Fo — L 111 ~r X 112 "T" L 221 "+" - 1 222 

-plm -plm _i_ pim _i_ pZm _i_ p/m ^99^ 

1 E ~ l 111 ~r L 122 "T" - 1 212 "+" L 221" l zz J 



The first three are the retarded three-point functions; T 1 ^ is retarded with respect to the 
first leg, T 1 ^ is retarded with respect to the second leg, and T 1 ^ is retarded with respect 
to the third leg. The other four vertices are related to the retarded ones through the KMS 



conditions J2TJ22 



r^ = JVi(r£ m - r 1 ™) + iv 3 (r^ m - r^ 
r& = N 2 (vi? - r£) + N 3 (rir - rS 

1 Fo - ^H 1 Ro - L R ) + ^2{i- R - 1 I 

T l ™ = r^f + T^ + T* R l ™ + N 2 N 3 (T% + Yf ) 

+ iv^r^ + r^ m ) + ivw^r^ + r^ m ) , (23) 



Ri) 1 
R ) 1 

Ri) > 



The vertex can be written as the outer product of three column vectors: 




(24) 



B. Integral equation for viscosity 

In order to perform a resummation we look at an integral equation of the form represented 
in Fig. [3] . In these diagrams, the vertex with the solid box is a full vertex whose form will 
be specified later. As we will see, it is the form of this vertex that determines which diagrams 
are being resummed. The four diagrams with corrected vertices in Fig. [3] are necessary in 
order to obtain a result for a that is pure real. Using corrected vertices on the right and left 
sides alternately is necessary to correctly take cuts of the vertex itself. Notice that diagrams 
with corrected vertices include an extra factor of r 3 which carries the index that corresponds 
to the momentum Q. The reason is as follows. The places where lines join in the bare one 
loop diagrams are not vertices and do not carry any factors, numerical or otherwise (expect 
for the kinematic factor i/ m ). In the diagrams with corrected vertices, the corrected vertex 



is given by (20). The factor of r 3 which is included in the definition of the vertex must be 
compensated for by an extra r 3 which carries the index that corresponds to the momentum 
Q. This extra factor of r 3 ensures that the diagram with the full vertex would reduce to the 
bare one loop diagram in the tree limit, which means that the combinatorics are correct in 
the integral equation that derives from Fig. [3] . The integral equation corresponding to the 
diagrams in Fig. [3a] is of the form, 

A^ = 2 j dKI lm (k,k)D ah (K + Q)D ha {K)I lm {k,k) (25) 

Using the decomposition (H) we obtain, 



(a _ a ') (1) = (A 21 - A 12 )W (26) 

= / dKIi m (k, k)[f k +q(a k - r k ) - f k {a k+q - r k+q )]I lm (k, k) 

where we define a k = DA(K),r k = Dr(K), f k = Dp(K) etc. Using the KMS condition ([?[) 
we obtain, 

(a - a') (1) = - J dK I lm (k,k)(N k - N k+q )(a k r k+q + h.c.)I lm (k,k) (27) 

In obtaining this result we have used the fact that we will eventually take the limit Q — > 



to obtain the viscosity fll2|). This limit gives rise to what is known as the pinch effect: terms 
with a product of factors a k r k+q and r k a k+q contain an extra factor ~ 1/Q relative to terms 
with products of propagators a k a k+q or r k r k+q . Thus terms proportional to a k a k+q or r k r k+q 
can be dropped. (The large terms occur when the contour is "pinched" between the poles 
of the two propagators, which gives rise to a factor in the denominator that is proportional 
to the imaginary part of these propagators) . 

The diagrams in Fig. [3b] correspond to an expression of the form, 

2Ai?= J dK t«T 1 c Z(K, Q, -K - Q)D db (K + Q)D bc (K)I lm (k J k) 

+ jdK I lm (k, k)D ac {K + Q)D da (K)T l ™(-K, -Q, K + Q)r b 3 (28) 

Using the decompositions (|S|) and ( |2~3D and the KMS condition (|7|) we obtain, 

(a-a') (2) = (A 21 -A 12 )( 2 ) (29) 

= \fdK{ a k r k+q [(N k+q - N k )T l ™ + N k+q T l ™ - N k T l ™ - Tf]I lm {k, k) 
- r k a k+q I lm (k, k)[(N k - N^T 1 ^ + N.T 1 ^ - N k+q T 1 ^ + T 1 ^} } 

where T x := T X (K, Q, —K — Q) and T^ := T x (—K, —Q, K + Q) and we have again used 
the pinch limit. We can further simplify this result by using the KMS condition (p3|) with 
the relations r^ (_) = [Tjg]* and r^ ( ~ } = -[T l p x ]*. We obtain, 

(a - a') (2) = \J dK {T l ™(K, Q, -K - Q)a k r k+q + h.c) (N k+q - N k )I lm (k, k) (30) 

Combining (|27D and (|30|) we obtain, 

a - a' (31) 

= JdK Uli m (k, k) + ^T%(K, Q, -K - Q))a k r k+q + /i.e.} (N k+q - N k )I lm (k, k) 

Notice that this expression depends only on the vertex T 1 ^ 1 . 

IV. VISCOSITY FROM CORRECTED FOUR POINT FUNCTIONS 

It is also possible to obtain a ladder resummation by writing the viscosity as an integral 
involving a corrected four-point function. This approach has advantages over the three-point 
function because it can be used to generate self consistent resummations. This point will be 
explained in detail below. 

9 



A. Notation 
The connected four-point function is given by the contour ordered expectation value, 

MC cd (X,Y,Z,W) = (t c mx)My)Mz)MW)) 

The 1PI four-point function is obtained by truncating external legs and forms a 16 compo- 
nent tensor which we can write as the outer product of four two component vectors, 



M 



: : ; 



The retarded 1PI four-point functions are given by 

M m = Mini + M m2 + M 1121 + M 1211 + M 1122 + M m2 + M 1221 + M 1222 
= - (x — y) (u + v) (w + z) (s + t) 



M1122 + M 



2112 



M2121 + M 



2122 



M K2 = Mim + M m2 + M1121 + M 2m 
1 
= -(x + y)(-u - u)(iu + z)(s + t) 

M m = Mim + M m2 + M 2m + Mi 2 ii + M 2112 + M m2 + M 2211 + M 2212 

= -(x + y)(u + v)(w — z)(s + t) 
M m = Mim + Mam + Mn2i + Mi 2 ii + M2121 + M2211 + M1221 + M 2221 

= -(x + y)(u + v)(tu + -z)(s - t) 

where we have used the relation 



(32) 



E M, 



abed 





a,b,c,d=l 




The other combinations are, 






M A = 


-(x + y)(u + v)(w — z) 


(s-t) 


M B = 


—(x — y)(u + v)(w — z) 


(s + t) 


M c = 


1 

—{x + y){u — v)(w — z) 


(s + 1) 


M D = 


-{x + y)(u — v)(w + z) 


(s-t) 


M E = 


—(x — y)(u — v)(w + z) 


(s + t) 


M F = 


1 

- (x — y) (u + v)(w + z) 


(s - 1) 


M a = 


— (x + y)(u — v)(w — z) 


(s-t) 



(33) 



(34) 
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M P = -j{x - y){u + v)(w - z)(s -t) 

1 
M 7 = —{x — y)(u — v )(w + z)(s — t) 

M s = -{x -y)(u- v)(w - z)(s + t) 

m t = -j[x -y){u- v)(w - z)(s- t) 

We use the decomposition of the four-point vertex: 

*('.)(:)C)(i)--C)('.)(i)('.)-"('.)C)(i)('.) 



B. Integral equation for viscosity 

The integral equation for the viscosity is obtained from Fig [4]. The one loop diagrams 
give, as before (|27|) . The second diagrams are of the form, 

X ab = / dP dK Ii m (p,p)Ii m (k,k) 

D ac (P + Q)D db (K + Q)D bf (K)D ed (P)M cdef (P + Q,-K- Q, -P, K) 

(36) 



We want to extract the combination A21 — A12 = (er — cr')^ 2 ). The result i 



is. 



1 f 

(a - a') (2) = - / dPdKI lm (p,p)I lm (k, k) [M R1 (r p+q f p (a k f k+q f k r k+q ) + f k f k+q r p+q r p ) 

—M R 2(f k a k+q (a p+q f p + f p+q r p )) + M R3 (f p+q a p (a k f k+q + f k r k+q )) 

+M m (f p+q r k (f p r k+q - f k+q r p ) - r k f k+q a p+q f p ) 

+M B (r p+q a p (a k f k+q + f k r k+q )) - M D (r k a k+q (r p f p+q + / p o p+g ))] (37) 

where we have suppressed the momentum arguments of the four-point functions. 

To simplify this result, we take the pinch limit and use the KMS conditions. These 
conditions for the four-point functions are analogous to the KMS conditions for the three- 



point functions (g^) and are derived in p2 |. The only one that we need is 



11 



(N k - N k+q )[M B - N P M R1 + N p+q M m ] = {N p - N p+q )[M* D + N k M* m - N k+q M* m ] 
Substituting fl38f) into Q37|) we get, 

^ {2) =-ljdPdKI lm (p,p)I lm (k,k)(N k -N k+q ) 



(38) 



[a — a 



[(M B - N P M R1 + N p+q M m )r k+q a k r p+q a p + h.c.) 



(39) 



Making the definition 



M{P + Q,-K- Q, -P, K) = M B (P + Q,-K- Q, -P, K) 

- N p M R1 (P + Q,-K- Q, -P, K) + N p+q M R3 (P + Q,-K- Q, -P, K) 

we can rewrite to obtain, 

1 /" 

(a - cr') (2) = -- / dPdKIi m (p,p)I lm (k, k)(N k - N k+q ) \Mr k+q a k r p+q a p + h.c. 

Combining (|27D and (|41|) we obtain from Fig. [4], 

a -a' = -J dK(N k - N k+q )I lm (k, k) { a k r k+q 

+ h.c. } 



(40) 



(41) 



Iim(k,k) + - J dPIi m (p,p)Ma p r p+li 



(42) 



V. LADDER RESUMMATIONS 

An argument similar to the pinch argument described in the section III-B gives rise to 
the conclusion that ladder graphs of the form shown in Fig. [5] contribute to the viscosity 
to the same order in perturbation theory as the bare one loop graph, and thus need to be 
resummed [jT0| , pT2|| . This effect occurs for the following reason. It appears that the ladder 
graphs are suppressed relative to the one loop graph by extra powers of the coupling, which 
come from the extra vertex factors that one obtains when one adds rungs (vertical lines). 
However, these extra factors of the coupling are compensated for by a kinematical factor 
~ 1/Q from the pinch effect, which becomes large when Q becomes small. One can also 
argue that the ladder graphs are larger than other graphs of the same order in the loop 
expansion. This conclusion is based on the following argument. We compare the graphs 
shown in Fig. [6a] and [6b]. The ladder graph shown in Fig. [6a] contains three pairs of 
propagators of the form D x D x+q which will give pinch contributions in the limit that Q — > 0. 
The graph in Fig. [6b] contains only two such pairs; the third pair is replaced by a pair of 
the form D x+y D x+q . The momentum variable y is an internal momentum which must be 
integrated over and thus there is a region of phase space for which it will be small. However, 
one expects that the measure will also be small in this region of phase space, and therefore 
that the contribution of all non-ladder graphs will be suppressed. This argument will be 
discussed in more detail in the next section. In this section, we will show that we can resum 
the ladder graphs by solving an integral equation for the corrected vertex. 
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A. Ladders from three-point vertices 

The ladder graphs are resummed by solving an integral equation for the three-point 
vertex. We can obtain T^ 1 self-consistently from the Schwinger-Dyson equation that corre- 
sponds to the diagrams in Fig. [7] . This expression for T 1 ^ can then be substituted into the 
integral equation (pT|) from which we obtain the viscosity. The diagrams in Fig. [7] give, 



r% c (K,Q,-K-Q) 

= J dPdRr b 3 r^I lm (p } p)D ca (R)D ac (K + R- P)D bc (P + Q)D ab (P) (43) 

+~ J dPdRT l ™ e (P } Q, -P - Q)D ad (P)D ec (P + Q)D ac (R + K- P)r^D ca (R)r^ 

We expand this expression by using the decompositions (§) and (|24j ) and performing the 
contractions. In the limit Q — > we make use of the pinch limit to obtain, 

T l ™(K,Q,-K - Q) = -^ J dP dR(I lm (p,p) + \r l ™(P,Q,-P - Q)) (44) 

[f p r p + q (f'r r + a'f r ) + a p r p+q (r'a r + a'r r + f'f r ) + a p f p+q (f'a r + r'f r )] 

where we have defined r' = Dr(P'); P' = R + K — P, etc. The key point is that this integral 
equation is decoupled: there are no vertex components other than T 1 ^ 1 on the right hand 
side. We introduce the notation, 

S(r, k)= a r r k + r r a k + f r f k 
0i (r, k)= f r a k + r r f k (45) 



(/>2(r,k)= f r r k + a r f k 



which allows us to write, 



T l ™(K,Q, -K-Q) = -^JdPdR(I lm (p,p) + \l%(P,Q, ~P - Q)) (46) 

[fpTp+gMr, r + k-p) + a p r p+q 6(r, r + k - p) + a p f p+q (j)2(r, r + k - p)] 



In conclusion then, we can obtain a resummation of ladder graphs by solving ( |4*6| ) for T 1 ^ 



l in 



and substituting into the expression (|31]) from which we obtain the viscosity. 



1. The Infrared Divergence 

It is easy to see that when the limit Q — > is taken, the resummation of ladder graphs is 
divergent. In fact, this result can be seen immediately from (|3i~D: the difference a — a' is of 
order Q° as Q — >• and therefore the viscosity, as obtained from ( |T^) and fll8|) , is divergent. 
This problem is customarily remedied by replacing the bare propagators on the rails (the 



horizontal lines) with hard thermal loop (HTL) effective propagators [l^,2^,p4j. We write 
the HTL propagators as 

D*(P) = - l - - (47) 
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where £(P) is the HTL contribution to the self energy. Note that the superscript * indicates 
a corrected propagator (and not complex conjugation). It is easy to see that this replacement 
regulates the divergence, but does not invalidate the conclusion that all of the ladder graphs 
need to be resummed. When resummed propagators are used, products of propagators of 
the form r* k+q a* k and r* k a* k+q contribute factors on the order of one over the imaginary part 
of the HTL self energy (instead of the factors ~ 1/Q that occur for bare propagators). The 
imaginary part of the HTL self energy is of order A 2 , and such a factor arises for each pair of 
rails. These factors will compensate the extra factors of A 2 in the numerator that accompany 
the addition of each rung. Since each additional rung is accompanied by an additional pair 
of rails, the result is that all of the ladder graphs have a piece that contributes to leading 
order. 

To resum ladder diagrams with HTL effective propagators, we rewrite fl31|) with HTL 
propagators, 



a — a 



(48) 
= -JdK(N k - N k+q ) [(I lm (k, k) + \v l ™{K, Q, -K - Q))a* k r* k+q + h.c.]I lm (k, k) 

and replace the vertex with the solution to the integral equation (compare (f46|)), 

r l ™(K,Q,-K-Q) = -^JdPdR(I lm (p,p) + ^r l g(P,Q,-P-Q)) (49) 

[/p*r; +(? 0*(r, r + k-p) + a* p r* +q S*(r, r + k - p) + a* p f* +q <p* 2 {r, r + k - p)} 

where 0*, 02, and S* are obtained from (^) with bare propagators replaced by HTL propa- 
gators. When the solution to ([4^) is substituted into ( f48"D we find that o — a 1 is of order Q, 
and the viscosity, as obtained from (ff2|) and (|18|) is finite. 

2. The Boltzmann Equation 

If we use the HTL self energy to two loop order and include a thermal mass and a 
thermal width, these results agree with results obtained previously in the imaginary time 



formalism |12| . In order to sum all planar ladder diagrams, Jeon has constructed a linear 
integral equation involving a 4 x 4 matrix valued kernel for the effective vertex, which is a 
four component column vector. The four components represent four different "cuts" in the 
imaginary time formalism. In the limit Q — > 0, the 4x4 matrix valued kernel reduces to 
a rank one matrix and the original matrix integral equation reduces to a decoupled single 
component equation for the reduced vertex, which is a complicated linear combination of 
the four different cuts of the vertex. We will show below that this result is considerably 
more straightforward to derive in the real time formalism. 
We define a new vertex 

In terms of this definition (H8f) becomes, 



o--a'~-JdK(N k - N k+q ) [T l ™(K, Q, -K - Q)a* k r* k+q + h.c.]I lm (k, k) (50) 
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The integral equation for Ti m is obtained by adding I\ m to both sides of one half of ( f49|) to 
obtain, 

r l ™(K,Q,-K-Q) = I lm (k,k)-^JdPdRT l ™(P,Q,-P-Q) (51) 

[f;r; +q {f>'* + kd + «;r; +g «r'* + r r v* + /;/'*) + a;r p+q (f:r>* + a* r n\ 



where we have used the definitions ([451) . 

We define the spectral density function as, 

p p = i(r* p - a* p ) (52) 

In the limit that Q is small the pinching terms dominate and we can rewrite the second 
term in fl5"T|) as, 



2 



->LjdPdRT%(P,Q,-P-Q) 

aK^a: + a'*r* r + f* r f*) - N p {f* r a'* + </'*) + iW/ r V* + </'*)] (53) 

Using (|47|) and (|52|) we have, 



and f p = —iN p p p . Taking Q — > and using these results we find that (|53|) becomes, 

- — / rfP di? —^——ppPrPp' [n r n p , +n r + n p (n r - ry )] (55) 

Using the identity 

n r n p i + n r + n p (n r — ry) = (1 + n p )(l + n p i)n r /{\ + n fc ) (56) 

and substituting into (|5i~|) we have, 

r^(^) = ^ - l^fc 2 -^-JdPdRdP' 

(2tt) 4 S\P + P'-R- K)p p p r pp, *\> {l + n p ){\ + n p ,)n r /{\ + n k ) (57) 

The appearance of the imaginary part of the HTL self energy in the denominator of this 
result indicates that a finite width is necessary to regulate an infrared divergence in the 
integral. Once the integral has been expressed in a regulated form, the imaginary part 
of the self energy can be dropped in the expressions for the spectral densities, since it is 
subleading. In (^) we take, 

p p -+ p° p = 2vrsgn(p )5(P 2 - rr&) (58) 

where mf h := m 2 + Re£#. We define 
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T l ™:=I lm (k,k)T% ear (59) 

hmV^i k)\ = T^-llm{k, K) = (klk m — — Ol m ) 
-n shear / r/-\ 

B{K) - toE^Tf)' 
and symmetrize to obtain from fl5"7|), 

T l g{K) = I lm (k, k)-—f dPdRdP' (2tt) V(P + P'-R- K) P y rP ° pl 

-(Ii m (p)B p + Iim{p')Bp> - Ii m (r)B r )(l + rip)(l + n pl )n r /(l + n k ) 
Rearranging we obtain, 

I lm (k, k) = ImZ R (K)B(K)i lm (k } k) + — -A , / dPdRdP' (2iv) 4 5(P + P'-R-K) 

12(1 + rifc) J 

p>° P ;,[i /m (p,p) J B(P) + I lm {p',p')B{P') - I lm {r,r)B{R)]{l + n p ){\ + n p ,)n r (60) 
In 4 theory ImE# can be expressed as, 

ImS R (K) = -— ( t — - — r ) / dP dR dP' (2vr) 4 (5 4 (P + P'-R-K) 
12 \(1 + n k )J J 

p pP y p/ (l + n p )(l + n pl )n r (61) 

which allows us to rewrite (|60|) as, 



I(k, k)im = *' > IdPdRdP' (2k)%P + P'-R- K)p pP y p , (62) 

12(1 + n k ) J y y 

[I lm {p,p)B{P) + I lm {p',p')B{P') - I lm {k, k)B{K) - I lm {r,r)B{R)]{l + n p )(l + n p ,)n r 



Equation (|62|) is an integral equation whose solution represents a non-perturbative re- 
summation of contributions to B. The last step, is to obtain an expression for the viscosity 



in terms of the quantity B. Using (|12l) , (|18D, ([50]) and (|59|) we have, 

rj = |- lim hm(l + n q ) f dK(N k - N k+q )k 2 [Tf ear (K, Q, -K - Q)a* k r* k+q + h.c.} (63) 

Using 

lim (1 + n q )(N k - N k+q ) = 2n k {l + n k ) (64) 

we obtain, 

77 = ^JdKk 2 n k (l + n k )[Tf ear alr* k + h.c.} (65) 

Using (JS3J) and ( |59"D and making the substitution p p — > p p we have, 
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n = ^JdKk 2 n k (l + n k )B k p° k (66) 

It is straightforward to show that our integral equation for the resummation of ladder 
contributions to the three point vertex (^) has exactly the same structure as the equation 
obtained from the linearized Boltzmann equation, and that our expression for the shear 
viscosity (|66|) is exactly the same as that obtained by defining the shear viscosity as a trans- 
port coefficient and using the linearized Boltzmann equation to get an integral equation for 
viscosity We begin from the Boltzmann equation for the phase space distribution func- 
tion f (X, P) which describes the evolution of the phase space probability density for the 
fundamental particles comprising a fluid, 

—^f(X,K) = \f dT 12 ^ k [fJ 2 (l + / a )(l + f k ) - (1 + /0(1 + f 2 )hf k ] (67) 

UJ k OX^ I J123 

where fi := f(X, KA, f k := f(X,K_) and dTi 2 ^ 3k is the differential transition rate for 
particles of momentum Pi and P 2 to scatter into momenta P% and K and is given by 

dT 12 ^ k := ^-\T{K,P z] P 2 ,P^W l= tg (2rr) 4 6(P l + P 2 - P 3 - K) (68) 

T is a multiparticle scattering amplitude describing a process in which particles of momenta 
Pi and P 2 scatter into momenta P% and K. Underlined momentum variables are on shell, 
and for the purposes of this section we specialize further to the positive mass shell. The 
Boltzmann equation is valid for distribution functions that describe the distributions of 
positive energy particles. 

In a weak coupling 4 theory, one can linearize the Boltzmann equation by expanding the 
non-equilibrium distribution function f(X,P) around a local equilibrium function fo(X,P): 

f(X,P) = f (X,P){l-x(X,P)[l + f (X,P)]} (69) 

where 

/o(X,£)=n(K5,)|; n(u) = — Ly (70) 



Expanding in powers of Vw one finds [12 



X(X,P) = P(X)A(X,P)V ■ u{X) + (3(X)B(X,P)[P ■ V(u(X) ■ P) - ^V • u(X)} (71) 

The coefficient A multiplying the divergence of the flow is related to the bulk viscosity. We 
obtain the shear viscosity from the coefficient B which multiples the shear in the flow. B 
satisfies the linear inhomogeneous integral equation fll0|| : 



x [l+n(wi)] [l+n(w 2 )] n(uj 3 )/[l+n(uj k )} (72) 

x (I lm (k, k)B(K) + I lm (p 3 ,P 3 )B(Ps) - I lm (P2,P2)B(P 2 ) - I lm { V i,Pi)B{P^) , 
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We want to compare this result with fl6"2"|). The two equations have exactly the same form 
once the delta functions have been used to do the frequency integrals in (p2|) and the classical 
piece, which corresponds to the positive mass shell, is extracted. The shear viscosity can be 
written in terms of B by looking at first order corrections to the energy momentum tensor, 
and comparing with the constitutive relation, 

(T V ) * -J^VjlVi(T°) + VjtfZ?) - j^VW>] - ^L^VW) + S^V) (73) 

where e := T 00 is the energy density, V = |T/ is the pressure, and 77 and ( are the shear and 
bulk viscosities respectively The result is [ H|r2 | 



ft r H 3 k 
ri = j~J 7^r^W(l + n(u k ))B(K) (74) 

15 J (2TT)' i UJk 



We want to compare this expression with our result (|66|). After doing the integral over the 
frequency components in ( |66|) and extracting the classical piece, the two expressions agree. 

B. Ladders from four-point vertices 

We can also resum ladder contributions to the viscosity by starting from the integral 
equation for viscosity in terms of the four-point function ( |3"§D and replacing the four-point 
function by the solution of the integral equation shown in Fig. [8]. The zeroth order 
terms come from the first graph on the right hand side of the figure. We obtain the term 



corresponding to the second graph in the figure by using the decomposition fl35l) . We have, 
A 2 



A r 

b = —r dR5(r, r + k — p) 

A 2 r 

— j dRdS [M B S{r, r + s- p)a s r s+q + M m (8(r, r + s- p)a s f s+q 



A 



M, 



+ 4> 2 {r,r + s - p)a s a s+q ) + M R1 (S(r, r + s - p)f s r 8+q + <f> x (r, r + s - p)r s r s+q ) ] 

A 2 r 
ri = -j- dR(j)i(r,r + k - p) 

A 2 r 
~ 32 7 dRdS [ M m($( r , r + s - p)r s r s+q + 0i(r,r + s - p)f s r s+q ) 

+ M m (f)i(r, r + s- p)a s f s+q + M B <f>i{r, r + s- p)a s r s+q ) ] (75) 

A 2 r 
M R3 = — j dR<j> 2 {r, r + k-p) 

A 2 r 
- — J dRdS [ M R3 (5(r,r + s - p)a s+q a s + 4> 2 (r,r + s - p)f s+q a s ) 

+ M m 4> 2 (r, r + s- p)r s+q f s + M B 4> 2 (r, r + s- p)r s+q a s } 
We can rewrite this result by taking the pinch limit and using the definition (fiDD. We obtain, 

M(P + Q, -K - Q, -P, K) = M(°» (P + Q,-K-Q,-P,K)-^JdRdS a s r s+q (76) 
M{S + Q,-K- S, -P, K)[5{r, r + s - p) - N p <j} X (r, r + s - p) + N p+q <f) 2 {r, r + s - p)] 



Note that as with ( fig ) this integral equation is decoupled. As a consequence, we can obtain 
a resummation of ladder graphs by solving (|76|) for M and substituting into the integral 
equation (pE2| ) from which we obtain the viscosity. 

VI. ANOTHER WAY TO RESUM GRAPHS 

In this section we will show that there is another way to obtain a resummation of con- 
tributions to the viscosity. We will see that the method discussed in this section indicates 
that the ladder graphs are not the only graphs that need to be resummed, in spite of the 
arguments discussed in the previous section. We will identify the extra graphs that are 
included in our resummation, which are not part of the conventional ladder resummation. 
The idea is as follows. In [2(J it is shown that in 3 theory, the retarded three-point vertex 



can be related to the self energy by a relation of the form, 

T R (K, Q, -K-Q)^ Q2 + l 2Q - - (U A (K) - U R (Q + K)) (77) 

where the symbol ~ indicates that we are looking at the infrared region of the momentum 
integrals in the various terms in the loop expansion of the three-point vertex and the self 
energy. This condition will be made explicit below. We will obtain an equation of this 
form for the corrected three-point vertices involved in (|31|), and for the corrected four-point 
vertices in (^). In both cases, a resummation is obtained by replacing the self energy 
in these expressions by the solution of a separate coupled integral equation that resums an 
infinite series of contributions to the self energy. This procedure is presented in detail below. 



A. The corrected three-point function 

In ([H]) the vertex (|20| ) is similar to the three-point vertex in <j) 3 theory except that it 
contains the composite operator iri m (X) (|i"4|). We will show that in the limit of infrared 
loop momenta this three-point function is related to a function that is similar to the sunset 
self energy from 4 theory. The equation that relates these quantities is similar to ([H]). At 
one loop, the relevant diagrams are shown in Fig [9]. In the vertex diagram, the place where 
the two lines join is a bare vertex T abc m fllTf). In the self energy diagram, a line carrying 



momentum P with an asterix contains a factor Ii m (p,p). Note that this notation does not 
indicate an insertion, but merely a multiplicative factor. We will show explicitly that the 
diagrams in Fig [9] satisfy the relation, 

Tf(K, Q, -K-Q)*> Q2 + 2 2Q - - (nf(K) - l4 m (Q + K)) (78) 

Note that in order to obtain this result we have rewritten Ii m (p+k,p + k + q) = Ii m (p + k,p + 
k) + Ii m {p+k, q) and taken only the first term. We have not included the extra contributions 
that result from the second term, since these terms will not contribute to the final result 
when the limit Q — » is taken. We verify this relation explicitly at one loop, and we also 
discuss the relation at higher orders in the loop expansion, working at zero temperature for 
simplicity, and discuss what sets of graphs are included. Finally, we will show that when 
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this expression for the vertex is substituted into the integral equation fl3~T|) from which we 
obtain the viscosity, and the self energy U lm is replaced by the solution of a separate integral 
equation, the result for the viscosity corresponds to an infinite resummation of graphs. 

1. The splitting relation for the three-point function 

We begin by deriving Q78|). The three-point function in Fig. [9a] is given by the expression 

r& = / dP J dR(-iX) 2 iD de (R)rliD' ed riiD eb (P)iD bd (P + Q)T b 3 I lm (p, P ) (79) 

We expand this equation using the vector column representation of the propagators (§J) and 
extract T 1 ^ 1 ( |22| ) . We make the change of variable P —>■ P + K. We will need the functions 



5(r, r — p), 4>i (r, r—p) and ^{r^r—p). These expressions are defined in (|45|). In the equations 



that follow we will suppress the momentum arguments for these quantities. We obtain, 
T l ™(K,Q,-K-Q) = -^JdPdRI lm (p + k,p + k) 

{fp+k r p+k+q4>l + fp+k+qO>p+k4>2 + l^p+k+q^p+k^) (80) 

Splitting the propagator pair and looking at the region of phase space where 

2P-Q^Q 2 + 2K-Q (81) 

we have, 



which gives, 

T%(K, Q, -K-Q) = - ^ - j dPdRI lm (P +k,p + k) 



t(g 2 + 2Q ■ k) 

[(flp+fe — r p+ k+q)S + fp+k<Pl — fp+k+q4>2] (83) 



We compare this result with the sunset self energy as shown in Fig. [9b]. We obtain, 

Il (K) l % = i l -j dRdP(t\) 2 i 3 T*r b 3 D ba (R)D ab (P + K)I lm (p + k,p + k)D ab (R - P) (84) 

Contracting indices and using the definitions of retarded and advanced self-energies (|Tl~D, we 
obtain, 



T1 (K)% = — - J dRdPI lm (p + k,p + k)(r p+k 5(r, r - p) + / p+fc 2 (r, r - p)) (85) 

A 2 r 
n (K) 1 ^ = -— J dRdPI lm (p + k,p + k){a p+k S{r, r-p) + f P+ k<f>i{r, r - p)) 

Comparing (|83D with ( 85|) gives the result (ff^). 
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When we look at (ffg) at higher order in the loop expansion, we find that the sets of 
graphs that satisfy this relation do not correspond to only the ladder graphs. Working at 
zero temperature for simplicity, it is straightforward to show that the vertex graphs shown 
in Fig. [10] correspond to the second order sunset graph shown in Fig. [11]. The important 
point is that there is a cancellation between contributions from the ladder graph (Fig. 
[10a]) and contributions from the two non-ladder graphs (Fig [10b, 10c]). This cancellation 
indicates that the leading contributions from the ladder graphs must be of the same order 
as the leading contributions from the non-ladder graphs, in contradiction to the prediction 
of the pinch argument, as described in the previous section. As explained previously, the 
conclusion that non-ladders are suppressed relative to ladders is based on the following 
argument. There is only a limited region of phase space (the infrared region) in which non- 
ladders are significant, and since the measure that corresponds to this region is small, the 
non-ladders should be suppressed. However, our derivation of ( f78|) is based on the realization 
that in the limit that Q — > it is precisely this region of phase space that is important, 
since both ladders and non-ladders are largest in this region. Since we have singled out the 
infrared part of the phase space, by imposing the condition 2P ■ Q <C Q 2 + IK ■ Q, the 
argument about the measure being small is no longer relevant and there are contributions 
from the non-ladder graphs that are as large as contributions from the ladder graphs. 



2. An Integral equation for II 

The second step in this procedure is to obtain an integral equation that resums an infinite 
set of contributions to II. The integral equation that resums the sunset contributions to the 
self energy is shown in Fig. [12]. We obtain, 

_A J dRdPT 3 a T*D ba (R)D ab (R - P)D ac (K + P)Tl l ™(K + P)D db (K + P) (86) 
Contracting indices we obtain, for the retarded and advanced combinations, 

H(K) 1 ™ = n (K) 1 ™ -^JdRdP (5[U(K + P)fr k+P r k+P ] 

+ 2 [n(iT + P) l ™f k+P r k+P + U(K + P) l ^a k+P f k+P + U(K + P)fa k+P r k+P }) 

U(K)f = U (K) 1 ^ - ¥■ J dRdP (S[U(K + P)fa k+P a k+P ] (87) 

+ 0![n(AT + P)ff k+P r k+P + U(K + P) l ^a k+P f k+P + T1(K + P)fa k+P r k+ 



where S, <f>±, 02 are functions f(r,r — p) and the momentum dependence has been sup- 
pressed. 

One of the advantages of the resummation technique outlined in this section is that it 
can be easily extended to include other types of graphs. By using an integral equation 
that resums different contributions to the self energy, one can obtain the resummation of a 
different set of contributions to the viscosity. 
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3. Contributions to the Viscosity 

Finally, we ask what contributions to the viscosity are obtained by substituting the 
solution of the integral equation (|87]) into (|78|) (which is in turn substituted into fl3T|)). The 
first term in (|3l| ) is the one loop graph. We will show that the second term is the first ladder 
contribution to the viscosity (Fig. [13a]), when we use (|78|) with the self energy given by the 
first order sunset graph (|85|) . Substituting (|78f) and fl35l) into fl3lf ) we obtain for the second 
term, 

X 2 f 1 /' 

(a-a'Y 2) = -— J dK -——— J dRdPI lm (p + k,p + k){(N k - N k+q )a k r k+q 

[{a P +k ~ r P+k + q )S + f p+k <pi - fp+k+qfo] + h.c.}Ii m (k, k) (88) 

As before, 5, (pi, <p2 are functions f(r,r — p) and the momentum dependence has been 
suppressed. We use fl8"2|) to obtain, 

A 2 
- / dK dtidf 

+<fiifp+kr p +k+q + faap+kfp+k+q) + h.c.}Ii m (k, k) (89) 



[a-a') {2) = ~y I dK dRdP I lm {p + k,p + k){(N k - N k+q )a k r k+q (5a p+k r p+k+q 



We want to compare this result with the result obtained from the first ladder graph shown 
in Fig. [13a]. Defining the quantity in the figure as X a b we want to look at a — a' = A 2 i — Ai 2 . 
From the figure we have, 

\ ab = A 2 J dRdPr c ^D cd {R)D dc {R - P)I lm {p + k,p + k) 

D da (P + K)D ac (P + K + Q)D cb (K + Q)D M (K)I lm (k, k) 

We substitute in the vector forms for the propagators @. We keep only the terms that give 
us non-zero contributions to (a — a')( 2 \ The result is, 

A 2 r 

(a - ct') (2) = — ;r / dKdRdPI lm (p + k,p + k){(N k - N k+q )a k r k+q (5a p+k r p+k+q + 

4>if P +kr p+k+q + 4>2a P +kfp+k+q) + h.c.}Ii m (k, k) (90) 

where we have made the change of variable K — > —K — P — Q in the second term in order to 
obtain a form that is easily recognizable as the hermitian conjugate. Once again, 5, <pi <p2 
are functions f(r,r — p) and the momentum dependence has been suppressed. Comparison 
of ( |8"9"D and ( PU[ ) verifies our statement that the first order sunset graph corresponds to the 
first order ladder graph. 

At higher orders the calculation is more complicated. In Appendix A, we discuss the two 
loop contributions, working at zero temperature for simplicity. We show that the second or- 
der calculation gives contributions to the viscosity of the form shown in Fig. [13b] and [13c] . 
The key point is that non-ladder graphs are included. The conclusion is that we can obtain 
a resummation of ladder and non-ladder graphs by substituting the solution of the integral 
equation that resums sunset graphs fl8~T|) into the expression for the vertex (|77[), which is in 
turn substituted into the integral equation fl3T|) from which we obtain the viscosity. 
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4- The Infrared Divergence 

The procedure described above will result in an expression that is divergent in the limit 
Q — > 0, in exactly the same way that the previous result, obtained by solving the integral 
equation for the vertex ([16|), was divergent. In this case, we expect to encounter difficulties 
in the limit Q — > since the condition (|81|) is not satisfied in this limit. We can remedy 
the problem by using HTL propagators. The vertex is the same as the vertex shown in Fig. 
[9a], except that the propagators are now HTL propagators fl4"T|). From ( |80| ) we obtain, 
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T l ™(K,Q,K-Q) = -^JdPdRI lm (p + k,p + k) 

\Jp+k r p+k+q ( Pl + Jp+k+q a p+k ( f ) 2 + r p+k+q a P +k" ) (91) 

We split propagator pairs as before. With HTL propagators, the splitting relation has the 
form, 

d* (p + k + Q)m (p + k)~ d-Ap + k)-d*r(p + k + Q) m 

In this case, the notation ~ indicates that we are looking at the region of phase space where 

2Q-P + P{Y>' A {K) - E' R {K + Q))<^Q 2 + 2Q-K- Y> R {K + Q) + E A (K) (93) 

This condition can be compared to the condition for the splitting of bare propagators (jBip . 



Since the imaginary part of the HTL sunset self energy is non-zero, the condition (|93|) is 
satisfied in the infrared region of the momentum integral, even when Q is taken to zero. 
The resulting expression has the form, 

A 2 / 1 \ 

r%{K, Q, -K-Q) = -— y Q 2 + 2Q . K + TiA {K)-Y, R {K + Q)) 

J dP dR I lm (p + k,p + k) [(a* p+k - r; +k+q )5* + f; +k <f>t - t p+k+q ^ (94) 

We compare this result with the sunset self energy (compare (|85|) ) . 

MK) 1 ^ = ~f dRdPI im(p + k,P + k)(r* p+k S* + /; +fe 0*) (95) 

MK)f = ~ J dRdPI lm (p + k,p + k)(a* p+k 5* + r p+k <P\) 

where S, <f>i, (fi 2 are functions f(r,r — p) and the momentum dependence has been sup- 
pressed. Comparing (0) and (|95|) we obtain, 

VtVCQ, -K - Q) = Q1 + 2Q . K _ E 2 R{Q + K) + SA(K) V'nK) - nS-W + K)) (96) 



We obtain the self energy from the integral equation (compare (|87|)), 
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IL(K)% = U (K)f -^JdRdP (5*[U(K + P)£V* +p r* +p ] 

+ ^ 2 [U(K + P)tf k+p rl + v + U ( K + P)>l +P r k+P + n(if + P)>* k+P r* k+P . 
IL(K)f = n (K)f-^JdRdP (8*[U(K + P) l Tal +p al +p ] (97) 

+ <PI[T1(K + P)tft +P rl +P + Il(K + P) l Tal +p r k+P + n(# + P)>* k+p r* k+p . 

The corresponding contributions to the viscosity are the same as in Fig. [13], but with all 
lines taken as HTL propagators. Note that, as usual when calculating with HTL propagators, 
we must include counter terms to avoid double counting. 



B. The corrected four-point function 

It is also possible to obtain an expression similar to (|77]) for the four-point function. We 
must look at the special case where two of the four fields have the same position co-ordinate, 
and the same Keldysh index. A vertex of this form is shown in Fig. [14]. Note that this 
restriction means that the 1PI four-point vertex has essentially the same form as a three-point 
vertex. We will show that the use of this vertex in the expression for the viscosity (|39| ) leads 
to a resummation of contributions to the viscosity, all of which contain one chain link. A 
diagram of this form is shown in Fig. [15]. It has been shown that diagrams containing chain 
links are suppressed fLOfl . The argument is as follows. As discussed in section [V-A-l] in the 



case of ladder diagrams, we expect that each additional pair of rails will introduce a factor 
1/A 2 from the imaginary part of the HTL self energy. Since each additional chain vertex 
introduces a factor A, it appears naively that chain diagrams contain a factor (1/A) n where n 
is the number of chain links. In fact this conclusion is invalid. A pair of rails in a chain link 
does not contribute a factor 1/A 2 because of the fact that the discontinuity of a chain link 
vanishes in the limit of zero external four momentum, and the real part of a bubblechain link 
does not contain pinching pole contributions. Consequently, chain diagrams are suppressed 
by a factor A n where n is the number of links. The four-point vertex is introduced only as 
a means to illustrate the self-consistent cancellation, as discussed in section [VII]. 

1. Integral equation for viscosity 

We need to express the viscosity in terms of an integral equation of the form shown in 
Fig. [16]. For the first term we obtain, 



A (1) 

A ab 



-i\ J dKdPI lm (k, k)I lm (k, k)D ac (P + Q)D ca (P)r^D cb (K + Q)D bc (K) (98) 

Performing the summations and taking the combination A21 — A^ 2 = {°~ — c') we obtain, 

(a - a') (1) (99) 

A /" 

= ~i- / dK dPIi m {k, k)Ii m {k, k)[(N p - N p+q )(N k - N k+q )(r p a p+q r k a k+q + h.c.)} 
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For the second diagram in Fig. [16] we obtain, 

J dK dP I lm (p, P)M(P + Q,-K-Q,-P, K) eced 
D ea (P)D ae (P + Q)D bd (K)D cb (K + Q)I lm (k, k) 



x(2) 
A ab 



Taking the appropriate combination we get, 

{a-o'f £ ) = \ fdKdPI lm (p,p) 
8 J 

{ (N k+q - N k )r k+q a k [M B (r p a p+q + a p r p+q + f p f p+q ) + {M m + M m )(a p f p+q + f p r p+q )] 
+ (N p+q - N p )r p a p+q (M B (f k f k+q + a k r k+q ) + (M D + M T )r k a k+q + 
(M m + M p )r k f k+q + (M R2 + M s )f k a k+q } I lm {k, k) (100) 



This result can be written in terms of a three-point function. Using (p3|), (|35D and (^) we 
find, 

M T = M D := iZ F] M B := 0; M m = M m := iZ R] 
M m = M 5 := iZ Ri] M m = Mfs := iZ Ro 

where the momentum arguments for the M's and the Z's are M X (P + Q, —K — Q, —P, K) 
and Z X (—K — Q,Q,K) and the vertices Z obey the same KMS conditions ( p3|) as the vertices 
Ti m . Rewriting the result in terms of the three-point function Z we find that the viscosity 
can be expressed in terms of one retarded three-point function: 

(a-a') (2) = l -JdKdPI lm ( P ,p)(N k -N k+q )(N p -N p+q ) (101) 

[a k r k+q a p r p+q Z R (-K - Q, Q, K) + h.c.}Ii m (k, k) 



Combining (PU|) and (|101|) we obtain 



a — a = — 
4 



JdKdPI lm ( P}P )(N k - N k+q )(N p - N p+q ) (102) 

(a k r k+q a p r p+q [-\ + Z R (K, Q, -K - Q)\ + h.c.)I lm {p,p) 



2. The splitting relation for the four-point function 

We consider the one loop diagrams shown in Fig. [17a,b]. The vertex shown in Fig. [17a] 
is given by, 

M bcha (P + Q,-K-Q,-P,K) 

= i?L J dSdRD bc (R + K + Q)D ab (R + K)D ac (S)D ca (R + S)r a r b T c 

Using the form of the propagators given in (Q) and extracting the vertex M R1 = iZ R we 
obtain, 
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A r 

Zr(K, Q, -K -Q) = — I dSdR [r r+k+q a r+k S + a r+k fr+k+q<Pi + f r+k r k+r+q (p2] (103) 

where 5, <p±, 02 are functions of the form f(s,r + s). Using (82) to split pairs of propagators 



of the form D r+k D r+k+q we obtain 

A^ l 

4 Q 2 + 2Q • K 



A f 

R = T r)2 i on . /^ / dSdR{a r+k 5 + f r +k4>2 - r k+r+q 8 - fk+r+ q 4>i\ 



where we are looking at the region of phase space where 2Q ■ R -C Q 2 + 1Q ■ K. We compare 
this expression with the results for the self energy shown in Fig. [17b]. We find that, 

Z « = - Qi + 2Q.K ^ (K) -^ K + Q)) (104) 

The lowest order contribution to the viscosity is shown in Fig. [17c]. Notice that the one 
loop graphs shown in Fig. [18] and the two loop graphs shown in Fig. [19] do not contribute 



to the expression ( 104 ) since the tadpole self energies are momentum independent. At higher 
loop order, the diagrams become increasingly complicated. Working at zero temperature for 
simplicity, the same relation is satisfied by the diagrams in Figs. [20-24]. 

3. An Integral equation for II 

In order to perform a resummation of a series of contributions to the viscosity, we need to 
obtain II as the solution of an integral equation which resums contributions to the self energy. 
As in section [VI- A- 2] , we consider the integral equation which resums sunset diagrams (see 
Fig. [13]). In this case, none of the lines contain crosses, because we do not need to work with 
a self energy that contains factors Ii m (p,p). As explained in detail below, this simplification 
is the motivation for introducing the four-point vertex. We find (compare ([87])), 

U(K) R = U (K) R -^JdRdP (5U(K + P) R r k+p r k+p 

+ 2 [n(K + P) R fk+ P rk+p + n(A" + P) A ak+pfk+p + U(K + P) F a k+p r k+p ]) 
U(K) A = U (K) A -^JdRdP (SU(K + P) A a k+P a k+p (105) 

+ 0i[n(if + P) R f k+p r k+p + U(K + P) A a k+p f k+p + IL(K + P) F ak+ P rk+ P }) 



4- Contributions to the viscosity 

As in section [VI-A-3], the last step is to ask what contributions to the viscosity are 
obtained by substituting the solution of the integral equation (|105| ) into ( |104| ) (which is in 



turn substituted into (|102|) ). It is straightforward to show that the diagrams produced are 
the same as for the case of the three-point vertex (as shown in Fig. [13]), except for the 
fact that each diagram contains one bare chain link insertion. One example is shown in 
Fig. [15]. From power counting arguments, it is known that diagrams containing chain link 
contributions are suppressed. We reiterate, that the four-point vertex is discussed only as a 
means to illustrate the self consistent cancellation, as discussed in the next section. 
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5. The Infrared Divergence 

As always, the result is infrared divergent unless the propagators are replaced with HTL 
propagators. At lowest order, the diagrams that satisfy 

Z R (K, Q, -K - Q) = - Q , + 2Q . K _^ R(Q + K)+ ^ A(K) ^K) - U R (Q + K)) (106) 

are the same as in Figs. [17a, b] with the propagators replaced by HTL propagators. Per- 
forming the substitution we obtain, at first order, the contribution to the viscosity shown in 
Fig. [25]. 



VII. SELF CONSISTENT RESUMMATIONS 

It has been shown that the resummation technique introduced in Section [VI] can be 
used to resum a larger class of diagrams than just the ladders. In the case of the corrected 
three-point vertex, the ladder graphs are included in the resummation, in addition to another 
large group of graphs that are shown to contribute at the same order. In the case of the 
corrected four-point vertex, we include the same large class of diagrams as before, expect 
that all contributions to the viscosity will have a chain link piece to the diagram. Now 
we want to consider a further increase in the set of diagrams to be resummed. We follow 
2"0fl . The basic idea is as follows. Until now, we have considered various forms of vertex 



corrections, with the propagator lines corrected by HTL self energies to obtain an infrared 
finite result. Now we will consider diagrams with corrected propagators of the form, 

6 < p > = p^m (107) 

where the self energies are not HTL self energies, but completely general expressions. Our 
goal is to find a cancellation between the self energies in expressions for the vertices ((|75|) 
and ( |104Q ), and the self energies in the corrected propagators. In [|(J, in a calculation of the 



self energy in 3 theory, it has been shown that these cancellations occur and are independent 
of the form of the self energy. The value of this approach is that the vertex corrections and 
the propagator corrections are written in the same form and dealt with simultaneously and, 
in this sense, they are treated self-consistently. 

We show below using the four-point function that it is possible to obtain a partial can- 
cellation of the type described above. We remind the reader that in addition to the fact 
that the cancellation is not complete, there are other problems with the four-point vertex 
approach. As explained in section [VLB], the decomposition of the four-point function in 
terms of self energies cannot be accomplished in general. We must restrict to four-point 
functions in which two of the fields are at the same space point, and the corresponding 
two Keldysh indices are equal. The consequence of this restriction is that we only obtain 
contributions to the viscosity in which there is one chain link piece. It seems likely that 
these two problems are related. We use the four-point vertex only as a means to illustrate 
the cancellation in the scalar theory. 

We begin with the integral equation shown in Fig. [26]. We have (compare ( |102j )), 
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a-a' = % -J dKdP(N p - N p+q ) {N k - N k+q )I lm (p, p) 

{a p+q r p a k+q f k (-X + Z R (-K - Q, Q, K)) + h.c.}I lm (k, k) (108) 

The propagators with bars are full propagators (|107|) , and the two bare propagators will give 



the bare chain link piece to each viscosity diagram, as discussed previously. Substituting 



in (|106|) we obtain, 



a — a 



% — JdKdP(N k - N k+q )(N p - N p+q )I lm {p,p)I lm (k,k) 



1 + Q1 + 2Q.K-Z r (Q + K) + Z a (K) {Ua{K) ~ Ur{Q + K)) ) °pW*^ + h - C - 

[Q 2 + 2Q-K- m R (Q + K) + 3U A {K) - Z R (Q + K) + Y, A {K)]a p r p+q a k f k+q + h.c 
Using D(P)- 1 = P 2 - n(P) we obtain, 

, y _ jb.f dK dP W - jW (^ - jW L(vv)L(kk) 

a a - t A J dK dP Q2 + 2Q K _ ^ {Q + r) + EA(K) a p r p+q l lm {p,p)l lm {k, k) 

(o fc - f k+q + [2U A {K) - 2U R {Q + K) + E A (K) - S fi (Q + K)]a k r k+q + h.c.) (109) 

When II is obtained as the solution of the integral equation (|105|) , Eqn. ( |109| ) represents 
the resummation of a huge number of contributions to the viscosity. Unfortunately, the 
cancellation between the vertex and propagator factors is only partial. We expect that the 
cancellation will be complete in a gauge theory with a three-point interaction like scalar 
QED. The motivation for this argument is discussed below. 

When calculating viscosity from the corrected three-point function, the presence of the 
composite operator (113) leads to difficulties. As shown in section [VI-A-1], this composite 



operator leads to a three-point vertex which cannot be decomposed in terms of a standard 
self energy, but rather the modified self energy which we have written as Hi m . Even after 
using translation invariance to write Il/ m = Iim.Il, there is in general no cancellation between 
the factor n from the vertex and the factor II from the full propagator. Note however that 
the numerical factor seems to be correct. Consider substituting (p6|) into (compare (pB), 



a-a' (110) 

= JdK {(I lm (k, k) + \v l ™(K, Q, -K - Q))d k f k+q + h.c.} (N k+q - N k )I lm (k, k) 

and replacing Ib m by Ij m II. If we could write II = fl then the cancellation between the 
factors of II in the propagators and the vertex would be complete: there would be no factors 
of II in the numerator of the resulting expression. 

We expect that the cancellation that we have described in this section does occur in a 
gauge theory like scalar QED as a result of the additional constraints imposed by gauge 
invariance. Note that equations of the form (|77|) look very much like Ward identities. It has 



been shown |2(| that the corresponding equation for scalar QED is 
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^=- ie wi^K^ A(K) -^ (K+Q)) (111) 

and corresponds to the statement that the part of the vertex that is large in the infrared limit 
must satisfy the usual Ward identity with the polarization tensor. In a future publication, 
we will look at this type of cancellation in viscosity in scalar QED. 



VIII. CONCLUSIONS 

We have introduced several different techniques for performing non-perturbative resum- 
mations. We used the closed time path to reformulate the traditional resummation of ladder 
graphs and obtained the same integral equation as that obtained previously in the imaginary 
time formalism. We have developed a technique which uses a pair of integral equations that 
gives resummations of both ladder and non-ladder graphs. We have argued that these non- 
ladder graphs contribute at leading order. Finally, we have shown how this resummation 
technique can be generalized to treat vertex and propagator corrections self-consistently in 
a way that gives rise to some cancellation between vertex and propagator corrections. We 
might expect that, in a gauge theory, the constraints imposed by the Ward identities will 
facilitate this cancellation. 

The various resummation schemes that we introduce in this paper become increasingly 
difficult to understand in terms of diagrams. It seems likely that a more physical understand- 
ing of these resummations can be obtained by interpreting the resulting integral equations 
in terms of an effective kinetic theory description. The equivalence of the integral equation 
corresponding to the ladder resummation and the Boltzmann equation was first discussed 
in |TI| . It should be possible to obtain a similar interpretation of the integral equations that 
correspond to the resummation schemes that we have introduced. 

ACKNOWLEDGMENTS 

This work was supported by the Natural Sciences and Engineering and Research Council 
of Canada (NSERC), and the National Natural Science Foundation of China (NSFC). 



29 



Appendix A 



We start from the integral equation involving the three-point function Ti m (|3"T|). We 
work at zero temperature and consider only the term that corresponds to corrections to the 
one loop graph. Note that, at zero temperature, it is not possible to regulate the pinch 
singularity using the imaginary part of the HTL self energy. We use the zero temperature 
case for illustrative purposes only, to simplify the discussion at higher loop orders. At zero 
temperature we have, 

a = 2jdK T lm (K, Q, -K - Q)iD(K)iD(K + Q) I lm (k, k) (112) 

We use the zero temperature version of (|78|), 

T lm (K, Q, -K-Q)~ Q2 + 2 - K (U lm (K) - U lm (Q + K)) (113) 

We make the definitions, 

a = K-P 1 -P 2 

a' = K + Q-P 1 -P 2 

b = K-P 1 -P 2 + P 3 + P 4 

V = K + Q-P 1 -P 2 + P 3 + P 4 

and we write D(a) = D a , etc. We make throughout the approximation 

Iim(a, a) ~ Iim(a', a') ~ Ii m (b, b) ~ Ii m (b', b') 

since there is no difficulty with taking quantities in the numerator to zero. 
We look at the self energy shown in Fig. [27]. We obtain, 

U lm (K) = ~ J dP 1 dP 2 dP 3 dP 4 D(P 1 )D(P 2 )D(P 3 )D(P 4 )D 2 a D b I lm (b,b) 
Substituting into (|112|) and ( |113| ) we have, 



A 4 r 1 

a = — \dKdP x dP 2 dP 3 dP 4 — 

2 J d Q 2 + 2( 



Q 2 + 2Q ■ K 
D(P 1 )D(P 2 )D(P 3 )D(P i )I lm (b, b)(D 2 a D b - D l2 D' b )D{K)D{K + Q) I lm {k, k) 

We can rewrite this equation by using the following identity 

D 2 a D b - D' a 2 D' b = D a D b (D a - D' a ) + D' a D' b (D a - D' a ) + D a D' a (D b - D' b ) 

The next step is to use the splitting relation for the propagators to write, 

q> + \q-k {d °- d '^ d ° d " etc ' 
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We obtain, 

\4 r 
a=—JdK dP 1 dP 2 dP 3 dP A D{P l )D{P 2 )D{P 3 )D{P 4 )D{K)D{K + Q) 

[I lm (a', a)D 2 a D' a D b + I lm (a' , a)D a D'*D' b + I lm (b\ b)D a D' a D b D' b ]I lm (k, k + q) 

which corresponds to the diagrams in Fig. [13b]. This procedure can also be extended to 
higher orders. It can be shown that the polarization tensor in Fig [28] corresponds to the 
contributions to the viscosity shown in Fig. 
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FIGURES 




FIG. 1. One-loop skeleton diagram for shear viscosity in cp 4 theory. 
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FIG. 2. Symmetric combination of one-loop diagrams for the shear viscosity in cj) 4 theory. 





(a) 





(b) 
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FIG. 3. Symmetric diagrams for the shear viscosity in terms of a corrected three-point function. 
The squre blob respresents the corrected three-point vertex. 
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FIG. 4. Symmetric diagrams for the shear viscosity in terms of a corrected four-point function. 
The rectangle respresents the corrected four-point vertex. 
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FIG. 5. The infinite series of planar-ladder diagrams in X(f> theory. 





FIG. 6. Five-loop ladder- (a) and non-ladder-graph (b) for the shear viscosity. 
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FIG. 7. Schwinger-Dyson equation for the corrected three-point vertex which includes ladder 
graphs in <p A theory. 
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FIG. 8. Schwinger-Dyson equation for corrected four-point vertex which includes ladder graphs 
in (j) 4 theory. 
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FIG. 9. (a) The three-point vertex including one ladder insertion; (b)the sunset self-energy 
that corresponds to (a); the asterix represents a factor I\ m . 
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( a) 



(c) 



FIG. 10. (a) The three-point vertex with two ladders; (b) (c) non-ladder diagrams which con- 
tribute at the same order. 




FIG. 11. The second order sunset graph that corresponds to Fig. [10]. 



FIG. 12. The integral equation that resums sunset contributions to the self-energy. 
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FIG. 13. (a) The first order and (b) the second order contributions to the shear viscosity. 




FIG. 14. Corrected four-point vertex with two joined legs. 
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FIG. 15. Typical ladder diagram for the shear viscosity with one chain. 



FIG. 16. Contribution one chain to the shear viscosity. Ihe black square represents the 
corrected four-point vertex. 





FIG. 17. (a) The four-point vertex and (b)its corresponding sunset diagram; (c)The one-chain 
and one-ladder diagram for the viscosity. 




FIG. 18. (a) The one-loop four-point vertex and (b) the corresponding self-energy. 
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FIG. 19. (a) Some two-loop four-point vertices and (b) their corresponding self-energy. 
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FIG. 20. (a) Three-loop four-point vertices with ladders and (b) the corresponding sunset 
diagram. 



FIG. 21. (a) Three- loop four-point vertices and (b) the corresponding self-energy. 





FIG. 22. (a) Four-loop four-point vertex and (b) the corresponding self-energy 
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FIG. 23. (a) Four-loop four-point vertices and (b) the corresponding self-energy. 
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FIG. 24. Four-loop four-point vertices and (b) the corresponding self-energy. 
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FIG. 25. The shear viscosity with one-chain link and one ladder. The box indicates an effective 
propagator. 
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FIG. 26. The shear viscosity with corrected 4-point vertex for diagrams with one chain link 
and effective propagators. 




FIG. 27. The second order sunset diagram. 




FIG. 28. The third order sunset diagram. 





FIG. 29. Some of the ladder and non-ladder diagrams corresponding to the third order sunset 
diagram. 
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